A two-sphere model for bacteria swimming near solid surfaces 
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We present a simple model for bacteria like Escherichia coli swimming near solid surfaces. It consists of two 
spheres of different radii connected by a dragless rod. The effect of the flagella is taken into account by imposing 
a force on the tail sphere and opposite torques exerted by the rod over the spheres. The hydrodynamic forces 
and torques on the spheres are computed by considering separately the interaction of a single sphere with the 
surface and with the flow produced by the other sphere. Numerically, we solve the linear system which contains 
the geometrical constraints and the force-free and torque-free conditions. The dynamics of this swimmer near a 
solid boundary is very rich, showing three different behaviors depending on the initial conditions: (1) swimming 
in circles in contact with the wall, (2) swimming in circles at a finite distance from the wall, and (3) swimming 
away from it. Furthermore, the order of magnitude of the radius of curvature for the circular motion is in the 
range 8-50 /jm, close to values observed experimentally. 

PACS numbers: 47.63.Gd, 47.63.mf, 47.55.dr 



I. INTRODUCTION 

In recent decades, interest in the dynamics of self-propelled 
organisms such as bacteria, fish, and birds has increased enor- 
mously. This can be explained, firstly, by interest in the 
physics behind the phenomena, for which tools from con- 
tinuum mechanics are used to describe the motion of sin- 
gle swimmers as well as from statistical physics to deduce 
collective behaviors due to their mutual interactions. Be- 
sides, study of how this locomotion affects biological pro- 
cesses such as human reproduction or bacterial infection re- 
mains fundamental 1 . Finally, development of applications in 
which artificial swimmers could perform specific tasks in mi- 
crofluidic devices is certainly an area of intense research 2 -. 

A particular kind of self-propelled organisms are bacte- 
ria. Considering their micrometer scale and typical propul- 
sion velocities, they are in the low-Reynolds-number regime, 
in which inertial effects are negligible in comparison with vis- 
cous ones. Therefore, the dynamics of the fluid is governed by 
the Stokes equation and immersed objects by the force-free 
and torque-free conditions 3 . In this category, the bacterium 
Escherichia coli (E. coli) has been intensively studied, and 
much is known about its genetics, biological processes, and 
motility 4 . The body cell is about 1 /im in width and 2-5 ftm 
in length, and around 70% of it is water. Its propulsion is due 
to rotary motors connected to around six filaments of approxi- 
mately 15 yL/m in length, which emerge from different points of 
the cell. When these filaments rotate counterclockwise, seen 
from behind, they form a bundle and the bacterium swims in 
an approximately straight line at about 20 /on/s-. This is called 
the run mode. If one or more filaments start to rotate in the 
opposite direction, the bundle (flagella) unravels, and the bac- 
terium can change its direction, which is known as a tumble. 
The alternation of these modes creates a three-dimensional 
random walk^. 

In addition, when E.coli swims near a non-slip surface, it 
describes clockwise circles when is viewed from above. In 
1995, Frymier et al£ reported a three-dimensional tracking of 
the movement of two kind of E.coli: wild type and a smooth- 



swimming mutants, which are unable to perform the tum- 
bling motion described above. In the first case the bacteria 
describe circular segments but they leave the surface after a 
tumbling. Conversely, the smooth-swimming cells swim for a 
much longer time along the surface. An explanation for this 
circular motion was given by Lauga et al. from the fact that 
the head and the flagella rotate in opposite directions, which 
together with a stronger drag near a wall, produce a net force 
directed to the centre of the circle. 

An extremely interesting question is how these living en- 
tities can modify the mass, momentum, and energy transport 
properties in a suspension. Since the pioneering work of Wu 
and Libchaber^, considerable efforts have been made to under- 
stand this issue. Recently, Mino et al.— observed that, near to a 
solid boundary, a suspension of E. coli increases the diffusion 
coefficient of passive tracers by a factor proportional to the 
active flux, i.e., the concentration of active bacteria multiplied 
by their mean velocity. 

In the literature, there are various models for E. coli with 
diverse levels of complexity; for example, Ramia et alJ2- sim- 
ulated E. coli as a sphere joined to a single helicoidal flag- 
ellum, and the velocity field near the surface was obtained 
using the boundary element method. In that work, the radius 
of curvature was about 10 /im, and it was reported a tendency 
for the bacterium to crash into the wall, i.e., it approaches 
the solid surface until it touches it, with the tail pointing up, at 
which point it stops. A different approach was made by Lauga 
et al~, who used the same bacterium geometry, but to obtain 
its trajectory near the surface, they calculated the resistance 
matrix for the swimmer. It is important to note that in their 
work the distance to the wall was fixed, not allowing vertical 
motion, which is an interesting point since the radius of cur- 
vature depends strongly on the distance to the wall. In 2008, 
it was reported by Li et alJ^, using fluorescence microscopy, 
that this gap is only tens of nanometers. Finally, it is impor- 
tant to mention the work of Gyrya et al.—, in which E. coli 
was modeled using two spheres joined to a point of force by 
a dragless rod (a similar model was introduced previously by 
Hernandez-Ortiz et al^). There, the velocity field was ob- 
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tained in the bulk for an isolated swimmer as well as for a 
pair of them. This latter approach inspired our work aiming 
to build a simple model for a bacterium using spheres, like in 
Gyrya et ai, but adding the interaction with the surface and 
rotation of the spheres, which is necessary to reproduce circu- 
lar motion near the surface. Construction of a minimal model 
of a bacterium that can perform this motion near a surface is 
the central aim of this work. 

The present work is organized as follows: in Sect. HI1 the 
point-force model for bacteria is described, and the corre- 
sponding problem close to a solid boundary is discussed. 
Then, in Sect. [Ill] our finite-size model is presented, obtain- 
ing its velocity in the bulk and the geometrical parameters of 
the model are fixed. The problem of this swimmer near a sur- 
face is treated in Sect. IIVI where the formalism of resistance 
tensors is introduced. After that, in Sect. [V] the results from 
numerical simulations are presented. Different types of mo- 
tion are obtained, depending on the initial conditions. Finally, 
conclusions are presented in Sect. [VT] 



II. BACTERIUM AS TWO POINT FORCES 
INTERACTING WITH A SOLID BOUNDARY 

At large distances, a pusher bacterium such as E. coli cre- 
ates a velocity field perturbation that can be described as 
that produced by a force dipole. Using this model, Berke 
et al*^ showed that a pusher swimmer experiences an attrac- 
tion toward surfaces in the form of a negative vertical velocity. 
Moreover, the swimmer orientation is affected by the pres- 
ence of the surface, adopting a stable equilibrium orientation 
parallel to the surface-^. However, the dipolar model fails to 
provide an appropriate description of the swimmer-surface in- 
teraction, as we show in the next paragraphs, due to the dom- 
inant finite-size effects from the swimmer. In particular, these 
effects regularize the apparently divergent vertical velocity ob- 
tained when the swimmer approaches the surface^. 

Using the model described in detail in Sect. [HI] in which the 
swimmer has a finite size, it is possible to compute the limit 
of point forces. Considering two point forces of intensities fo, 
separated at a distance L, analytic expressions can be obtained 
using the images systemic. If the distance from the swimmer 
to the surface is h (Fig. [TJ and its orientation is parallel to the 
surface, the vertical velocity is given by 



V 7 = - 



3/oL/i 3 



27T77(L 2 +4/2 2 ) 5 / 2 ' 



(1) 



where 77 is the dynamic viscosity of the fluid. 

For a dipole of vanishing size (L — » 0) and fixed force 
dipole of strength p = foL, the vertical velocity matches the 
result of Berke et ai, which diverges when h — > 0. Note that 
limits do not commute and the velocity divergence is only 
obtained when the limit in L is taken before the limit in h. 
However, when L is finite, the velocity is regularized for all 
distances and vanishing for h — > 0. Furthermore, the angu- 
lar velocity vanishes as predicted by Berke et aM. The factor 
(L 2 + Ah 2 ) 1 ! 2 corresponds to the cross distance between one 
sphere and the image of the other. Therefore, the attraction to 




FIG. 1. A two-point-force model for a bacterium near a wall. The in- 
teraction with the wall is obtained using the image method proposed 
by Blake and Chwangi^. The arrows represent the forces exerted by 
the point forces on the fluid for a pusher bacterium. Dashed lines 
represent the interactions due to the image forces. 



the wall is due to the two-sphere interaction through the sur- 
face. In fact, it can be easily verified that the field created by 
the image of a point force moving in the x direction does not 
have a vertical velocity or angular velocity at the position of 
the point force itself and therefore cannot lead to an attraction 
toward or repulsion from the surface. 

It must be remarked that the point particle approximation 
is not accurate, as it predicts a vanishing vertical velocity 
when the height of the force dipole vanishes (h — > 0), while 
the bacterium should stop its vertical motion when its body 
reaches the surface. Near-field corrections, obtained from 
the lubrication-theory approximation, should be used to de- 
scribe the swimmer motion accurately. The model described 
in Sect. [HI] considers this approximation, and in Sect. [V] the 
full swimming dynamics is described. 



III. THE TWO-SPHERE MODEL 



The bacterium is modeled as two spheres, denoted by H 
(head) and T (tail), with radii a# and aj, respectively, and 
connected by a dragless rod of length L (Fig.[2j. The positions 
of the spheres are r# and r^, and the unit vector along the 
swimmer is defined as n = (rH-rrf/^H-rrl = (rfj-r^/L. To 
mimic the effect of the rotary motors in E. coli, the rod applies 
equal and opposite torques +Tb = ±Ton on the spheres, with 
positive sign for H. The T-sphere represents the bacterium 
flagella, and in the present model its rotation is not directly 
related to the propulsion; rather it is considered as an extra 
propulsion force foh applied by the fluid at the sphere center. 
The rotation of the spheres along the swimmer axis and the 
interaction with the solid surface are the main differences from 
the dumbbell models proposed befor e) 12,13 . 

Due to the micrometer scale of the swimmer, the total force 
and torque acting on each sphere are zero. Making explicit 
the separation between forces and torques applied by the fluid 
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A. Swimmer in the bulk 
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FIG. 2. Two-sphere model for bacteria: two spheres of different radii 
connected by a dragless rod, which exerts equal and opposite torques 
±Tb- Over the T-sphere, a propulsion force /oft is applied by the fluid 
to mimic the effect of the rotating flagella. 



(hydro) and those applied by the rod, 





= F^ ydro + F int = 0, 


(2) 




= F !; ydro - F int = 0, 


(3) 




= ^+^ = 0, 


(4) 


t t 


= T h / dl "° -T B = 0, 


(5) 



where F lnt is the internal force exerted by the rod. The hy- 
drodynamic forces and torques are those over smooth spheres 
with the nonslip boundary condition (for example the Stokes 
drag in the case of spheres in the bulk). The exception is the 
/"-sphere where, to mimic the propulsion effect of the rotating 
flagella, an additional propulsion force /oft is added to F^ dl °. 
From Eqs. dU and (0, 



(t h - t t ) ■ ft = (r^ dro - r h / dro ) ■ n + 2t b ■ n = 0. (6) 



The torque intensity to = Tb ■ n has been measured for E. coli, 
obtaining values of the order of 1 pN jjm (to = 1 .26 pN /im in 
Ref value that we use in the present article). 

From the geometrical constraints of the swimmer, the an- 
gular velocities of the spheres, perpendicular to the swimmer 
axis, must be equal (denoted by £2#), and since the rod is rigid 
in these directions, 



In the bulk, neglecting the Faxen corrections 3 , the hydrody- 
namic forces and torques on the spheres are given by 
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These expressions consider the far-field approximation, that 
is, when the separation between the spheres is large in com- 
parison with the radius. The neglected contributions are of 
order 0(a 2 /L 2 ). Here u h (th) is the velocity of the fluid at the 
position of the center of sphere H but calculated in its absence 
(the same being valid for sphere T). Being more specific, if 
the swimmer is isolated, the field produced by the two spheres 
in the far-field approximation is 



u(r) = -G(r-r w ).F^ dro . 



G(r-r r )-F h r ydl "°, 



(13) 



where G(r) is the Oseen tensor defined as G(r) = 
+ ^p-) with I the unit tensor. Again, the Faxen cor- 
rections 0(a 2 /L 2 ) haven been neglected. Note that F hydr0 are 
forces applied by the fluid. Hence, the fields u H (r) and u r (r) 
are 

u w (r) = u(r) + G(r - r fl ) ■ F^ ydl "°, (14) 
u r (r) = u(r) + G(r - r T ) ■ F^ 5 " 110 . (15) 

To obtain the velocity of the swimmer in the bulk, like in 
Ref.— , we replace Eqs. (TBI and ( fl3T l in d9UT2l. and with 
the zero-force and zero-torque conditions together with the 
geometrical constraints, we obtain 



U„ = U 7 



2L - 3a H 



12m][L(ciH + aj) — 3fl//fl7-] 

„3 



a\-I? 



Ton, 



T n. 



/oft = U n, (16) 



(17) 



(18) 



(£l H -fi r )xn = 0. (7) 



Besides, the linear velocities of the spheres are related by 



U H = V T + Ln x Sl B , (8) 



where Uh (Ut) an d (&r) are the linear and angular veloc- 
ity of the sphere H (T). 



Thus, the linear velocity is proportional to the propul- 
sion force, with a constant of proportionality depending on 
the swimmer geometry. Similarly, the angular velocities of 
the spheres are along the swimmer axis (i.e., it swims in a 
straight line), with magnitude proportional to the imposed 
torque along the axis, depending also on the relation between 
the geometrical parameters an, aj, and L. Since the torques 
applied by the rod have opposite sign, the //-sphere rotates in 
clockwise direction while the /"-sphere rotates counterclock- 
wise, both as seen from behind the swimmer. 
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B. Choice of geometrical parameters 

The geometrical parameters of the model aj, and L 
were chosen in the following way: the head sphere was re- 
lated to the cell body using the equivalent sphere concept^. 
Following Lauga et al, an equivalent sphere having the same 
viscous resistance as a prolate ellipsoid moving along its axis 
of symmetry has radius 



a = 



4w( 2(f> 2 -l (> + V^ 2 - 1 ' 



(02 _ 1)3/2 



In 



1 



(19) 



where w is the width and <p is the aspect ratio of the ellipsoid. 
In the case of E. coli, the radius of the equivalent sphere varies 
from 0.89 to 1.07/im 7 . We therefore fix an — 1 yum. 

Recent experiments have shown that, when the flow field 
around an E. coli bacterium is measured, the separation of 
the associated force dipole is on the order of 2 yum^l. With 
this information, we fix the separation between the spheres at 
L — 2 fim. 

We took the radius of the tail sphere as aj = 0.5 pm. This 
choice is rather arbitrary, since this model is a great simplifi- 
cation of E. coli and a direct comparison is not possible. Other 
values give qualitatively similar results. 

Finally, equation ( TToT l allows one to fix fi, for a given ge- 
ometry by knowing the experimental value obtained in Refi^, 
Uo - 20/mi/s. Using the values above, fo = 1.13 pN. 



IV. TWO-SPHERE SWIMMER IN THE PRESENCE OF A 
SOLID BOUNDARY 

When the swimmer is near to a solid boundary, there is a 
complex hydrodynamic interaction with the wall. One possi- 
ble approach to calculate this interaction is to use the image 
method proposed by Blake and Chwangi^, i.e., to consider the 
bacterium as two equal and opposite point forces and point 
torques. Since this is a far-field approximation, and the swim- 
mer can be very close to (namely touching) the wall, we will 
introduce corrections given by the Faxen formula^ and lubri- 
cation theory-^. 

Using the resistance matrix formalism, which is based on 
the linearity of the Stokes equations, the hydrodynamic forces 
and torques on the spheres that compose the swimmer are 
given by the relation 
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where the 6x6 matrices Kh and %r give the interaction of 
each sphere with the wall, and 'Rut and Kth are the couplings 
between them considering the wall. These matrices already 
take into account the Faxen and lubrication corrections. In the 
article of Swan and BradyJ^, expressions for these matrices 
are given in the case of spheres of the same radius in the far- 
field approximation; the generalization to spheres of different 



radii is rather complicated. This approach, however, is too 
rigid, because if two or more swimmers are considered, the 
size of the resistance matrix grows accordingly, and even more 
so the complication of the computation of its elements. 

As a simplification to this approach, we compute the inter- 
action of the swimmer with the wall as the complete interac- 
tion of each sphere with the boundary (considering the far- 
and near-field contributions), while the coupling between the 
spheres is obtained from the interaction with the external flow 
produced by the other sphere. For example, for the //-sphere, 



F hydro 
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hydro 
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H 
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H 

(21) 

where the matrix elements % are the 3x3 resistance tensors, 
and (U^, £2^) correspond to the fluid velocity and fluid an- 
gular velocity at the position of the //-sphere due to the pres- 
ence of the /"-sphere. This approach has the benefit of being 
simple enough to extend to the case of many swimmers or 
other passive objects, or to include an imposed flow. The as- 
sociated velocity fields are introduced in the (U°°, £1°°) terms. 
This approximation is valid when the two spheres are suffi- 
ciently separated and no close field interactions are present. 
As a consequence, the corrections are of order 0(a 2 /L 2 ). 

The hydrodynamic coupling between the spheres is ob- 
tained by computing the flow produced by each of them and 
evaluating it at the position of the other sphere. In the presence 
of a plane wall, the Oseen tensor is modified by image terms. 
The image set of a force monopole is not simply an opposite 
force as in electrostatics, but also includes a force quadrupole 
and a mass source dipole; the combination of these together 
with the original force monopole produces the flow that satis- 
fies the nonslip boundary condition at the wall— In the case 
of a sphere close to the plane, the flow should be modified to 
include the Faxen correction of the spheres and the lubrication 
flows. In the present model we will neglect the Faxen correc- 
tions, as they contribute subdominant terms. Thus, the flow 
produced by the //-sphere (that applies a force F and a torque 
r over the fluid), evaluated at ri = (x\,yi,z\), is 
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(22) 



where y = \y\,y2,h} is the position of the force point, r = 
{xi - y\, X2 - yi, xt, — h] is the vector from the singularity to 
the observation point, and r = [x\ —y\,Xz— yi, x 3 + h) is the 
vector from the position of the image to the observation point. 
In the coordinate axes, x and y are the planar directions and z 
is the perpendicular direction to the plane, as shown in Fig.Q] 
Considering the symmetries and the linearity of the flow at 
low Reynolds number, the resistance tensors have the follow- 
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ing general expressions for a single sphere: 



ular to a wall in the far-field approximation i: 
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and similarly for the tensors that couple with the velocity field, 
K FU °° , <R FCr ° , <R tU °* , and <R Tn °* . 

In the bulk, the resistance tensors take the asymptotic values 
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where a is the radius of the corresponding sphere. Together 
with the geometrical constraints, these produce the bulk re- 
sults <[T6]), OS, and (US]). 

To obtain each term in the resistance tensors for a sphere in 
the presence of a planar wall, we take advantage of the linear- 
ity of the Stokes equation. In fact, it is possible to consider 
separately the cases of a sphere moving and rotating in a qui- 
escent fluid and that of a stationary sphere in an ambient fluid. 
The first case allows <R FU , H Fa , 7? Tt/ , and <R Til to be computed, 
and the second case is used to obtain ft" 7 " , <R m ° , 7? Tf/< \ and 



A. Translational and rotational motion of a sphere in the 
presence of a plane wall in a fluid at rest at infinity 

In the far-field limit, the resistance for a sphere of radius a 
translating parallel and perpendicular to a wall, with its center 
located at distance h from the wall (Fig. [3), is 17 
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(26) 



Furthermore, due to the presence of the solid boundary, there 
is a coupling between the translation and the torque 
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where 7? Fn = ['R tU ] t (see Refi^ for a description of the gen- 
eral properties of the resistance matrix). Similarly, the resis- 
tance experienced by a sphere rotating parallel and perpendic- 
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FIG. 3. Sketch of the quantities used to compute the resistance ten- 
sors. The surface is in the x-y plane, the sphere radius is a, h is the 
distance from the sphere center to the plate, and S = (h — a) is the 
gap- 



On the other hand, when the sphere is close to the plane 
[that is, when the gap distance is much smaller than the sphere 
size, 6 = (h- a) <sz a], the asymptotic lubrication theory must 
be used. Goldman et aL- obtained the following asymptotic 
expressions: 
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(29) 



We propose global interpolation formulas, valid in both 
limits, with the following generic form 



[fl + b x] log 



) + c + do- +e () — + f () —, (30) 



where x = h/a is the normalized distance from the center of 
the sphere. The coefficients of the interpolation are fixed to 
match the two asymptotic expressions up to the known orders 
in each case. For example, for 'R FU , the expression is 
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1 .5296 - 0.9963 -| log (^—^ 
a J \ h 



1.9963, 



+ 0.4689 g) + 0.4327 0.1358 (f) . (31) 

The expressions of the global fits in the other cases are given 
in AppendixlAl 

On the other hand, for a sphere moving perpendicular to the 
wall, we use the resistance obtained by Bevan and Prievei^ as 
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a regression of the exact solution given by Brenner— 
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(32) 



It is remarkable that this expression is valid for both short and 
long distances from the wall. Finally, rotation perpendicular 
to the wall has been rarely treated in the literature; in the work 
of Lee and Leal^I, the Stokes equation is solved for a sphere 
moving and translating in the presence of a wall. The expres- 
sions for the resistances are given in terms of infinite series. 
In particular, 
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where 770 = -arccosh(jt) and c = yjx 2 - 1, with x - h/a. 
For small gaps, this expression goes as fFf 1 = 8^77a 3 [ 1.202 + 
0.414 (x - 1) log(jc - 1)]. A global interpolation between this 
limit and the expression d28l ) is proposed in the form 
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and the coefficients are given in AppendixlAl 



B. Immobilized sphere in an ambient flow 

In the work of Goldman et al. that considers a shear flow 
passing around an immobilized sphered, it was found that the 
shear flow induces forces and torques. The associated resis- 
tances approach finite values when the sphere is in contact 
with the wall: 
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As in the case of a sphere in the bulk, for large distances, the 
resistances are the same as for a sphere translating and rotating 
in a quiescent fluid but with opposite sign, 
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(36) 

The global interpolation between both limits is given in Ap- 
pendix A. 

On the other hand, when the flux is perpendicular to the 
wall, there is no known expression for the lubrication approxi- 
mation. For simplicity, we assume that, as for ftf 1 u and 9C. u " , 
this resistance does not diverge when 6 — > 0. Therefore, we 
use the expression for large distances of a sphere moving per- 
pendicular in a rest fluid, with the opposite sign, i.e., 
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Finally, there remains the coupling between force and rota- 
tion (and equivalently torque and translation) due to this ex- 
ternal flow. In Ref. — these coupling were not considered, be- 
cause the imposed flow was a simple linear shear flow charac- 
terized by a single parameter, the shear rate. Hence the local 
velocity U°° is proportional to the local angular velocity 
We have extended this to consider external flows that may not 
be linear and that are characterized by the local velocity and 
angular velocity, independently. A condition for consistency 
is that, in the case of linear flows, the known result should be 
recovered. For simplicity, and without any other information, 
we neglect these couplings. 



V. RESULTS 
A. Stable swimming in contact with the surface 

At each time step, the linear system of equations for the 
linear and angular velocities is solved numerically. The asso- 
ciated matrix depends on the orientation and vertical position 
of the swimmer, becoming singular when any of the spheres 
touches the wall. The swimmer motion, that is its position and 
orientation as a function of time, is obtained using a fourth or- 
der Runge-Kutta method with adaptive time step. 

For a wide range of geometrical parameters (0.5 yum< an < 
2 fim, 0.2 yL/m< aj < 1 /jm, and 1 /jm< L < 4 /urn), includ- 
ing those in Sect. MI Bl the dynamics of our model is as fol- 
lows. If the swimmer is initially pointing toward the surface, 
it approaches the surface in almost a straight line until it is 
very close to it, at which point it slows down and orients such 
that both spheres almost touch the surface (close but not ex- 
actly parallel to the surface because of the difference in radii). 
Then, the complete swimmer is attracted to the surface until it 
touches it in finite time. At this stage, the swimmer stops any 
motion (translation and rotation) in a stable state except for 
eventual numerical errors if the time step is not small enough. 
This motion is similar to that previously reported for other 
models: when the initial conditions of the swimmer are pre- 
pared such that it reaches the surface, it crashes into the sur- 
face and stops moving 10 . In our model, the swimmer stops 
with both spheres touching the surface. 

The stopping process can be understood because the re- 
sistance tensors diverge when the spheres touch the surface. 
Therefore, the finite propulsion force and rod torque cannot 
produce motion. However, the nature and character of the di- 
vergence of the different tensors is different, requiring deeper 
analysis. The vertical resistance fl™ d32l diverges strongly 
as l/<5. This divergence is necessary to avoid penetration of 
the body into the surface and it is always present, quite in- 
dependently of the precise geometrical shape of the swim- 
mer. On the other hand, the parallel resistance QTT l has 
a soft divergence as log(<5), which is only noticeable when the 
sphere is extremely close to (namely, in contact with) the sur- 
face. However, before becoming divergent, in real swimmers, 
the roughness at the nanoscale regularizes this resistances 3 -. 
The same happens for all resistances with the exception of 
'R FU . Therefore the resistances are naturally regularized at the 
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nanometer scale, and eventually electrochemical forces enter 
into play, also avoiding that the spheres touch the surface^. 

To keep the model at a hydrodynamic level, without includ- 
ing extra forces, the vertical resistance < R F ^ U is modified such 
that it diverges at a small but finite value of 6; therefore, the 
vertical movement is stopped before the other resistances di- 
verge. Specifically, instead of fl32l i, we take 

K v 2 + 9(f) + 6(tf 

7 — i — > 

6*na 2 (f) + 6(f) 

where 6* = (6- e) and the gap size is chosen as e — O.OIah = 
10 nm. This value agrees with the swimming distance of some 
tens of nanometers measured by Li et al., which is consis- 
tent with the prediction of the stable position including the 
Derjaguin-Landau-Verwey-Overbeek (DLVO) potential 1 1 . 

With this regularization scheme, the swimmer can perform 
stable motion in contact with the surface, as described in the 
next section. 



B. Swimming modes for fixed parameters 

Depending on the initial condition, characterized by the ini- 
tial height of the sphere H, z(0), and the initial orientation with 
respect to the surface a(0) (Fig.|4^), the swimmer shows dif- 
ferent asymptotic behaviors. Three different behaviors were 
identified: (I) swimming in circles in contact with the plate, 

(II) swimming parallel to the plate at a finite distance, and 

(III) swimming away from the surface. Figure [4}5 presents a 
phase diagram of the initial conditions that lead to these final 
regimes, which are schematically presented in Fig. |4j;. Ani- 
mations of these three regimes can be seen in the Supplemen- 
tary Material^. 



1. Swimming in circles touching the surface 

When the swimmer is initially close to the surface and 
pointing towards it (negative angle a), it approaches the sur- 
face, orients such that both spheres almost touch the surface, 
and swims in circles. This state is stable, and the swimmer 
remains in this state almost in contact with the surface thanks 
to the regularization described in Sect. IV A| All initial con- 
ditions in region I of the phase diagram end in the same fi- 
nal state, characterized by zh - Rh - Zt - Rt - £■ Con- 
sequently, the final orientation, which depends on the swim- 
mer geometry, is given by a = 14°. The swimmer velocity is 
Vdrcie = 9.8 jum/s, and the radius of curvature Rc = 15.8 /urn. 
Experiments show a wide dispersion of velocities and ra- 
dius of curvature of E. coli swimming in contact with a solid 
surface^. The computed values for this model fall within the 
observed experimental ranges. 



2. Swimming in circles at a finite distance 

When the initial orientation angles are higher than in the 
former regime, while still not being too large, the swim- 
mer ends up in a stable state, moving parallel to the sur- 
face at a finite distance. The final configuration is given by 
Zh = 1-51 fim, zt = 0.84 yum, and orientation with respect 
to the surface a = 19°. Due to the distance to the wall, the 
hydrodynamic interaction with it is strong, and the radius of 
curvature of this state is Rc = 12.9 /urn while the swimming 
velocity is that of the bulk V pa raiiei = 20 jum/s. All the initial 
conditions in region II of the phase diagram end in the same 
state. This regime has not been identified in experiments as a 
different regime from the previous one. Possibly because the 
radii of curvature are similar and the distance to the wall is not 
easy to measure. 

3. Swimming away from the surface 

Finally, if the swimmer starts with orientations even larger 
than those of the previous regime, it will swim away from 
the surface. It must be remarked that, if z(0) is large, even if 
the initial angle is negative, the reorientation produced by the 
interaction with the wall can be strong enough to take it away 
(not shown in the phase diagram of Fig. |4j. The final states in 
this regime are not unique and depend on the initial condition. 

C. Sensitivity on the gap distance 

In Section [VA"l it was argued the need of a small gap e, that 
would take into account the roughness at the nanometer scale 
and its value was fixed to e = 10 nm. The motion in the bulk 
and the swimming mode parallel to the surface are insensitive 
to this parameter choice not like the circular motion in con- 
tact with the surface. It is found that when changing the gap 
in the range 4 nm< e < 20 nm, the radius of curvature and 
the circular velocity vary in the range 8 //m< Rc < 50 //m and 
8.6/iin/s< Vdrcie < 10.8 yum/s, respectively. This sensitivity on 
the value of e could be responsible of the wide dispersion of 
velocities and radius of curvature observed in experiments^, 
due to the diversity in the nanometric roughness of the bacte- 
ria. 



VI. CONCLUSIONS 

A simple model for a swimmer at zero Reynolds number 
has been built, consisting of two spheres of different radii 
joined by a dragless rod. The model aims to reproduce the 
swimming of a pusher bacteria like E. coli near solid surfaces. 
The hydrodynamic forces and torques on each sphere are com- 
puted taking into account the full interaction with the surface 
using the complete resistance matrix. The hydrodynamic in- 
teraction between the spheres is treated approximately by con- 
sidering the force and torque over a sphere under the flow pro- 
duced by the other. The rod imposes opposite torques on the 
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a) b) c) 




z(0) [urn] 



FIG. 4. (a) Initial condition for the swimmer: the orientational angle with respect to the surface is o-(O), and the position of the center of sphere 
H is z(0). The geometrical parameters are a H = 1 fan, a T = 0.5 fan, L = 2 fan, and the gas size is set to e = 0.01 fan . (b) Phase diagram of 
the initial conditions. Three situations are observed: (I) swimming in circles in contact with the wall, (II) swimming parallel to the wall at a 
finite distance, and (III) swimming away from the wall. The final states for regimes I and II are indicated by a symbol. The upper left corner 
is forbidden by the condition that the T sphere must be above the surface, (c) Sketch of the three final regimes. 



spheres making them rotate in opposite directions. Finally, 
the propulsion is obtained by adding an extra hydrodynamic 
force on one sphere (the tail) as to mimic the effect of the ro- 
tary flagella of E. coli. The dimensions of the spheres and 
their separation are fixed according to experimental results of 
the size of E. coli, its velocity in the bulk, and its effective 
force dipole. 

When pure hydrodynamic interactions of perfect spheres 
are used, if the swimmer is set to approach the surface, it 
reaches the surface in finite time with both spheres in con- 
tact with it. At this point the swimmer stops moving due to 
the divergence of the hydrodynamic resistances. Appealing 
to the existence of nanometric roughness of the body cell or 
electrochemical forces with the surfaces, the resistances are 
regularized at the nanometric scale obtaining finite velocities 
when it swims close to the surface. 

The swimmer shows three different behaviors depending 
on the initial condition. If the swimmer is close to the sur- 
face pointing toward it, it approaches the surface and orients 
in such a way as to have both spheres touching the surface. 



Once reached this state, the swimmer performs a circular mo- 
tion at constant speed with a radius of curvature in the range 
8-50 jum, depending on the value used to regularize the re- 
sistances, which is in the range of the experimental observa- 
tions. In the present model, the swimmer remains in this state 
forever. It is expected that surface roughness, thermal noise, 
velocity agitations or the tumbling mechanism can allow the 
swimmer to exit the wall region. A second regime is obtained 
when the swimmer starts at higher altitudes. In this case, it 
approaches the surface without touching it and swimming in 
circles parallel to the surface with a gap of approximately a 
half of micron. The trajectory is also circular with a similar 
radius of curvature but larger velocity. To our knowledge this 
regime has not been identified experimentally. Finally, if the 
swimmer is initially far enough, it will reorient as to move 
away from the surface. 
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Appendix A: Fits for the resistance tensors 



Using the notation shown in Fig. [3] the normalized distance to the wall is defined as x — h/a, the gap distance is 6 — (h—a), and 
the minimum distance that the swimmer can approach to the wall is e. Global fits are built such that the appropriate series when 
6 <K a and when x » 1 reproduce the known asymptotic expansions (see discussion in Sec. IIVb . There is not global bounding 
error of the fits, but we expect that the interpolating procedure produce uniformly convergent expressions if more known terms 
were added in both limits. The resulting resistances for a sphere moving in a quiescent fluid and for a sphere immobilized in an 
ambient fluid are the following. 
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1. Global fits for a sphere translating and rotating in a quiescent ambient fluid 
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2. Global fits for a sphere immobilized in an ambient fluid 
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